Radio Emission from an Ultraluminous X-Ray Source 
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The physical nature of ultraluminous x-ray sources is uncertain. Stellar mass black holes with 
beamed radiation and intermediate mass black holes with isotropic radiation are two plausible 
explanations. We discovered radio emission from an ultraluminous x-ray source in the dwarf 
irregular galaxy NGC 5408. The x-ray, radio and optical fluxes as well as the x-ray spectral 
shape are consistent with beamed relativistic jet emission from an accreting stellar black hole. If 
confirmed, this would suggest that the ultraluminous x-ray sources may be stellar-mass rather 
than intermediate mass black holes. However, interpretation of the source as a jet-producing 
intermediate-mass black hole cannot be ruled out at this time. 

One of the most enigmatic results to emerge from x-ray population studies of galaxies is the discovery of non- 
nuclear x-ray sources which appear to have luminosities factors of 10 to 100's times the Eddington luminosity 
for a neutron star - the so-called ultraluminous x-ray sources (000- The physical nature of these sources is 
controversial and possible models include very luminous supernova remnants (4), gamma-ray burst remnants 0, 
accreting stellar-mass compact objects with radiation which is beamed mechanically (6) or relativistic ally Li 51 . or 
accreting intermediate mass black holes (2i with masses between 100 to several 1000 solar masses. The detection 
of x-ray variability in several sources suggests that these sources are powered by accretion of material from a 
companion star or interstellar gas QISJ ar, d are not supernova or gamma-ray burst remnants. If the radiation from 
these variable sources is isotropic, then a high inferred luminosity requires an intermediate-mass black hole. If the 
radiation is beamed, then a stellar-mass black hole could explain the observations. 

Optical counterpart searches have shown that ultraluminous x-ray sources preferentially occur near star-forming 
regions ( 9 1 although they are also found in old globular clusters U0ill\ . In star-forming systems, one unique mas- 
sive star counterpart has been identified ( 72 ) and nebulae have been detected around several sources U3t . In some 
cases, counterpart searches have shown that the sources are background objects (74). Further identification of 
optical and radio counterparts is needed to reveal the physical nature of the ultraluminous x-ray sources. 

The dwarf irregular galaxy NGC 5408 contains a well studied ultraluminous x-ray source, 2E 1400.2-4108 
J75ll76ll7^ l2l. We observed the galaxy for a total exposure of 4683 s with the Chandra X-Ray Observatory 
on 7 May 2002 using the High-Resolution Mirror Assembly and the Advanced Camera Imaging Spectrometer 
spectroscopic array operated in an imaging mode suitable for a high count rate source il <St . Images in the 0.3- 
8 keV band showed a point source (we place an upper limit on the source diameter of 0.9") at a position of 
a(J2000) = 14h 03m 19s.63, <5(J2000) = -41° 22' 58".7 with an error radius of 1" at 90% confidence and a flux 
of 2.7 x 10~ 12 erg cm -2 s -1 in the 0.3-8 keV band corrected for Galactic absorption along the line of sight. If 
the x-ray emission is isotropic, then the source luminosity would be 1.1 x 10 40 ergs -1 in the 0.3-8 keV band at 
the distance of 4.8 Mpc to NGC 5408 The long-term x-ray light curve (Fig. 1) shows variability at a level of 
40%. 

A 4.8 GHz radio map of NGC 5408 from the Australia Telescope Compact Array (ATCA) showed a radio 
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source at the location of the x-ray source. We reprocessed the 8 hour ATCA observation of NGC 5408 taken on 3 1 
March 2000 (20 1 to produce a natural weighted map that provides optimum signal to noise ratio for point sources 
at the cost of angular resolution. We detected a 4.8 GHz radio source with a flux of 0.26 ± 0.04 mJy (la error) 
at a position of a(J2000) = 14h 03m 19s.606 and <5(J2000) = -41° 22' 59".572 (J2000) with an error radius 
of 0.1" at 90% confidence i21\ . consistent with the Chandra position within the astrometric uncertainties. No 
source is present at the same location in 8.64 GHz data obtained simultaneously with the 4.8 Ghz data and we 
place an upper limit (3a) on the flux of 0.12 mJy based on the rms noise near that location. From the 4.8 GHz 
detection and the 8.64 GHz upper limit, we place a lower bound on the radio spectral index a r > 1.0, defined so 
S v oc v~ ai '. From radio source counts i23i at 4.9 GHz, the probability of chance occurrence of a 0.18 mJy or 
brighter background radio source within the Chandra error circle is less than 1.0 x 10~ 5 . To estimate the probability 
of chance coincidence with a radio source internal to the galaxy, we compare with NGC 1569 which is a dwarf 
irregular similar to NGC 5408 and has two sources which would have a 4.8 GHz flux of 0.18 mJy or larger if 
placed at the distance to NGC 5408 (22 1. The probability of chance coincidence within 1" of the Chandra source 
given 2 radio sources within the optical extent of NGC 5408 is less than 0.002. Further, all of the radio sources in 
NGC 1569 have a r < 0.6, which is inconsistent with the steep spectrum of our radio source. 

Several optical sources in a Hubble Space Telescope (HST) image of NGC 5408 (Fig. 2) lie within the 
ATCA/HST error circle; note that the relative error circle is dominated by the HST astrometric uncertainty. The 
brightest object in the radio or x-ray error circles has a visual magnitude of 22. 1 . The colors of the objects detected 
in both bands are similar to those of the overall stellar population of the galaxy M 9\ . None of the objects show 
evidence of extended emission. 

Given the radio counterpart and the brightest optical counterpart in the Chandra error circle (which we take as 
an upper limit to the optical flux), we find that the x-ray to optical flux ratio \24\ is f x / f v > 380. We can also 
calculate two-point spectral indices 125II26> between the radio, optical, and x-ray bands using the flux density at 
4.8 GHz, the V band, and 1 keV. We find that the radio to x-ray index a rx = 0.16, the optical to x-ray index 
a ox < 0.18, and the radio to optical index a ro > 0.15. The ratio f x /fv of the object is much larger than the ratios 
found for stars, normal galaxies, clusters of galaxies, and AGN (24 26 I. We rule out the possibility that the object 
is a member of any of these classes of source. Identification as a blazar is ruled out because the f x J f v exceeds 
that of known blazars (26 27), the radio spectrum is inconsistent with those of blazars which have a r < 0.8 (27i, 
and the x-ray to radio flux ratio is 9 x 10~ 9 which is an order of magnitude larger than the most extreme known 
blazar 1271 . 

Some young supernova remnants have high x-ray luminosities. However, supernova remnants produce copious 
optical line emission, which is not detected here. Taking the brightest source in the Chandra error circle, the upper 
limit on the Ha emission is 7 x 10~ 15 ergcm~ 2 s _1 if all of the flux in the HST F606W filter is Ha emission. A 
similarly obtained upper limit on the [Olll] emission is 1.6 x 10~ 14 ergcm~ 2 s" 1 . These upper limits are below 
the line fluxes measured for highly x-ray luminous supernova remnants (28 29). Optical line emission could be 
hidden by a high obscuring column (301, but such obscuration is inconsistent with the relatively low absorbing 
column density found from the x-ray spectrum (Fig. 3). Also, the observed x-ray variability would be unexpected 
for a supernova remnant. The line fluxes and the f x / f v ratio are also inconsistent with identification of the source 
as a young supernova. 

X-ray binaries and rotation-powered neutron stars sometimes have high f x / f v ratios similar to what we find 
for the source in NGC 5408. The x-ray variability indicates the source is not a rotation-powered neutron star. If 
the source was located in the Milky Way, the luminosity would be of order 10 33 — 10 34 ergs -1 . Radio emissions 
when detected from Galactic x-ray binaries in such low luminosity states i31\\32\ have a flat or inverted spectra 
which is inconsistent with the steep spectrum of the source and suggests that the source is not a low-luminosity 
x-ray binary in the Milky Way. The steep radio spectrum is consistent with optically thin synchrotron emission as 
produced by x-ray binaries in high luminosity states or during state transitions i33\\34\ [32 1 and is consistent with 
the source being an x-ray binary in NGC 5408. The source x-ray variability may be less than seen from many x-ray 
binaries, but is consistent with that of Cyg X-l. 

Radio emissions from x-ray binaries are associated with relativistic jets. The radio emission we detected 
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from the ultraluminous x-ray source in NGC 5408 is probably from a relativistic jet. It could be from an accreting 
stellar-mass black hole producing relativistic jets (a microquasar) beamed toward us ( 33 I. The intensity of isotropic 
rest-frame emission from a continuous jet with a power law spectrum with radio spectral index a r is amplified by a 
factor S 2+ ° r where 6 = 7 _1 (1 — PcosO)^ 1 is the Doppler factor (for the approaching jet) where (3 is the jet speed 
as a fraction of the speed of light and 7 = (1 — /3 2 )~ - 5 . The Galactic microquasar GRS 1915+105 has 5 = 0.34 
and a distance near 11 kpc {34). GRS 1915+105 is highly variable. Radio fluxes of about 300 mJy at 4.8 GHz are 
observed in bright radio flares when superluminal jets are resolved, whereas the more typical radio fluxes are about 
10 mJy while the source is x-ray active but not undergoing a bright radio flare. If the ATCA observation happened 
to catch the NGC 5408 source during such a bright flare, then a Doppler factor as low as 8 = 2.1 would suffice to 
produce the observed radio flux (with a r = 1.0). However, the probability of catching such an event with a single 
radio observation is low unless the NGC 5408 source is an extremely active jet source. Scaling from the more 
typical radio flux of GRS 1915+105 of about 10 mJy, we find that a Doppler factor S > 5.8 would be needed to 
produce the observed flux. This is compatible with the constraints on the ejecta velocities for GRS 1915+105. If 
7 = 5 as inferred for GRS 1915+105, then tilts of the jet axis relative to the line of sight of up to 10° could produce 
S > 5.8. Roughly 1 out of 70 of extragalactic microquasars should be aligned like this. A steeper radio spectral 
index or higher intrinsic luminosity would reduce the required S. 

The x-ray emission could also be relativistically beamed jet emission l?5t . The fact that a r > a rx rules 
out the possibility that a single component spectrum, like synchrotron emission, produces both the radio and x-ray 
emission. Inverse-Comptonizationof photons from a massive companion star by a jet may power the ultraluminous 
x-ray sources \36\ . The x-ray spectrum predicted by this process has the form of a broken power law with the break 
near 1 keV and a cutoff at energies above the x-ray band. The Chandra x-ray spectrum is consistent with this form 
(Fig. 3) with a break energy of 0.65 ± 0.06 keV and a spectral index at high energies of a x = 2.1 ± 0.2. The 
beaming pattern for inverse-Compton emission with a photon source external to the jet, e.g. from the companion 
star, varies as S 3+2oia: for a continuous jet, so very large amplifications are possible, e.g. 4 x 10 5 for S = 6. Hence, 
the jet x-ray luminosity may be as low as 10 34 ergs" 1 . The break energy we measure is similar to what has been 
modeled i36b . However, our measured x-ray spectral parameters would imply an unusually steep power law index 
for the injected relativistic electrons of p = 4.2 ± 0.4. 

If the source is an accreting black hole and the x-ray emission is isotropic, then the luminosity implies a 
minimum mass of 80M Q . The spectral index of the hard spectral component, if the spectrum is fitted with the 
sum of a power law plus thermal accretion disk emission, is a x = 1.75 ± 0.34 which is similar to that found for 
black hole candidate x-ray binaries in their high flux states {37s. The temperature of the soft thermal component 
is kT = 0.11 ± 0.03 keV which is lower than found for stellar mass black holes in the Milky Way i38t . Such a 
low temperature is also unlikely to arise from mechanically beamed radiation from a stellar-mass black hole ©. 
Some other ultraluminous sources have high temperatures, which is inconsistent with the presence of intermediate 
mass black holes unless they are very rapidly rotating (7). The low temperature found for the NGC 5408 source 
is consistent with that expected for a non-rotating black hole with a mass above 100M© QE). The radio spectral 
index and the radio to x-ray flux ratio are similar to the (highly variable) ones measured for microquasars i33\\34\ 
suggesting that the radio emission could arise from jets similar to those produced by microquasars. The source 
lies about 12" from the star formation regions of NGC 5408 that contain super star clusters, giving a projected 
displacement of 280 pc. An x-ray binary with a speed of 10 km/s, from either interactions in a cluster or a 
supernova kick, could traverse this distance in 30 My. Hence, an intermediate mass black hole could have formed 
in one of the observed super star clusters i39t and moved to the present location within the life time of a massive 
companion star. 

The ultraluminous x-ray sources have been interpreted as evidence for intermediate mass black holes 0. The 
results presented here show that the radio, optical, and x-ray properties of the ultraluminous x-ray source in NGC 
5408 are consistent with beamed emission of a relativistic jet from a stellar-mass black hole. This suggests that, at 
least some of, the ultraluminous x-ray sources may be beamed emission from stellar-mass black holes (135b . 
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Figure 1: Long term x-ray light curve for the x-ray source in NGC 5408. Data shown as diamonds are, in 
chronological order, from the Einstein IPC \40\ . the ROS AT all-sky survey, the ROS AT HRI M7\ . two observations 
with the ROSAT PSPC 03, the ASCA SIS (0, and Chandra. The cross marks the time of the ATCA observation 
(the flux is arbitrary). The absorbed flux in the 0.5-2 keV band, calculated from each instrumental count rate using 
a spectral model derived from fitting the Chandra data, is plotted. The source appears variable with a ratio between 
the maximum and minimum fluxes of 1 .4. 
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Figure 2: An optical image of NGC 5408 showing the region containing the ultraluminous x-ray source with the 
ATCA 4.8 GHz radio contours (in red at flux densities of: 0.12, 0.16, 0.2, 0.3, 0.6, 1.2, and 2.4 mJy/beam), error 
circles for the Chandra (blue) and ATCA (green) sources, and the main star formation regions (marked in blue as 3 
and 4, note the associated radio emission (20); 1 and 2 are foreground stars U7Y ). The arrow points North and has 
a length of 2". To allow for the various astrometric uncertainties, we used a 1.5" radius circle around the Chandra 
position and a 0.8" radius circle around the ATCA position to identify possible optical counterparts. Two HST 
WPFC2 exposures (EB, both 600 s and obtained on 4 July 2000, were analyzed using the HSTphot (41 1 stellar 
photometry package to obtain simultaneous photometry in the two filters (F814W and F606W) and corrected for 
reddening using an extinction E(B — V) = 0.068 based on dust maps (42 1 and an Ry = 3.1 extinction curve. 
The image shown is with the F606W filter which includes the lines Ha and [OIII]. We corrected the absolute 
astrometry of the HST images using stars from the USNO A2.0 catalog and imwcs tool from the Smithsonian 
Astrophysical Observatory Telescope Data Center. There are no USNO A2.0 stars on the WF3 chips where the 
x-ray source is located, and we estimate a total astrometric error for positions on the WF3 chip of 0.5". 
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Figure 3: X-ray spectrum for the ultraluminous x-ray source in NGC 5408 in the 0.4-8 keV band. The x-ray 
spectrum is inconsistent with simple single-component models. A broken powerlaw model with an absorbing 
column density equal to the Galactic line-of-sight value of 5.6 x 10 20 cm -2 gave a good fit with x 2 /DoF = 
62.5/61, a low energy powerlaw photon index between -0.8 and +1.7, a break energy of 0.65 ± 0.06 keV, and a 
high energy powerlaw photon index of 3.1 ± 0.2. The spectrum is also adequately fitted (x 2 /DoF = 62.6/60) 
with a model consisting of thermal emission from an accretion disk at a temperature kT = 0.11 ± 0.03 keV with 
an added powerlaw component with photon index T = 2.75 ± 0.34 and an absorption column density Nn = 
(3.5 ±1.9) x 10 21 cm -2 , significantly above the Galactic column. Here and in the text, all error quoted on spectral 
parameters are 90% confidence for a single parameter. The spectral modeling was done with the ISIS software 
(http://space.mit.edu/ASC/ISISl to correct for the moderate (14%) pileup due to the high counting rate. 
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